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Quasiparticles in a Bose-Einstein condensate are sensitive to space-time distortions. Gravitational 
waves can induce transformations on the state of phonons that can be observed through quantum 
state discrimination techniques. We show that this method is highly robust to thermal noise and 
depletion. We derive a bound on the strain sensitivity that shows that the detection of waves in the 
kHz regime is not significantly affected by temperature in a wide range of parameters that are well 
within current experimental reach. 


The detection of gravitational waves [l| remains an open problem and represents one of the most ambitious en¬ 
terprises of science in the 21st century. After years of active efforts mm, several large-scale experiments are still in 
operation around the globe, based both in laser interferometry -such as advanced LIGO, GEO 600 or VIRGO- and 
Weber bar detectors -such as AURIGA and Mario Schoenberg. However, no successful observation of a gravitational 
wave has been reported yet. Therefore, together with new upgrades of the existing setups, new major international 
projects are expected to start operations in the short and medium term, including space-based laser interferometers 
such as DECIGO and LISA. This enormous investment of resources is backed up by indirect evidences of the existence 
of gravitational waves as well as a number of experiments confirming the predictions of Einstein’s General Relativity, 
a theory from which the existence of spacetime ripples is a natural consequence [4]. However, since the Earth is very 
far from typical sources of gravitational waves, the intensity of the latter is so tiny when it reaches our detectors that 
gravitational wave detection is always a daunting task. |38| 

Recently, a new way of detecting spacetime distortions was proposed using a different physical principle [5]. The 
state of the quasiparticles of a Bose-Einstein condensate (BEC) is modified by the passing of the gravitational wave. 
If the frequency of the gravitational wave matches the sum of the frequencies of two BEC modes the transformation 
of the state is resonantly enhanced in a phenomenon resembling the Dynamical Casimir Effect 00 , characterised 
by a linear growing in time of the transformed state. Indeed, the scheme is equivalent to an artificial modulation of 
the length of the BEC trap, which can be implemented by a modulation of the atomic interaction strength 0 El .This 
gravitational quantum resonance is absent in laser interferometers since the frequency of a gravitational wave is very 
far from the optical regime and is also different of the vibrational resonances in Weber bars. In [5], we showed that 
the sensitivity of the setup is low enough to, in principle, enable the detection of gravitational waves in a certain 
experimental parameter regime. 

In this work, we show that our scheme (see Fig. 0 ) for gravitational wave detection is highly robust against 
the effects of thermal noise. In [5] it was assumed that the quasiparticle state is prepared in a particular pure state 
before the passing of the gravitational wave. Now we include an initial temperature and show that the sensitivity 
is not significantly affected in a wide regime of temperatures well within experimental reach. Thus, we give a step 
further in the analysis of the experimental feasibility of our scheme for gravitational wave astronomy. Following the 
same spirit, we also show that the spacetime ripple does not generate additional thermal depletion on the atomic 
bulk of the condensate, but only induces a phase shift. Indeed, this phase shift is much less sensitive to the action of 
the gravitational wave. However, there have been proposals to detect gravitational waves with atom interferometer 
setups, both Earth and space-based mm ■ These experiments would aim at detecting low-frequency gravitational 
waves, from 10 -5 to 10 -2 Hz. Our setup operates in a completely different frequency range, ranging from Hz to KHz 
-the same frequency range as LIGO. Moreover, the sensitivity improves with the frequency while the sensitivity of 
LIGO is optimal in the range of 100 Hz and decreases at higher frequencies. In particular, this implies that models of 
spacetime waves generated by the merging of neutron star binary systems into massive neutron stars and black holes 
go beyond LIGO’s capabilities while being in principle within the reach of our proposal [12 . Detecting gravitational 
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FIG. 1: Sketch of the setup: A BEC in a box-like potential that acts as a cavity for the phononic excitations. Two phononic 
modes m and n are initially prepared in a thermal two-mode squeezed state represented by the initial covariance matrix ao. A 
gravitational wave of amplitude e transforms the state producing excitations. The new state a e depends on e. Measurements 
on the modes can be used to estimate the amplitude of the spacetime distortion. 


waves in the kHz regime would allow us to deepen our understanding of neutron stars by gathering information about 
their mass and radius. These parameters are necessary to describe the Neutron star state equation[13] and also allow 
cosmologists to compute distances that are key in the study of the cosmological constant and dark matter [IT. Other 
recent proposals of high-frequency gravitational-wave detectors can be found in [T5l [IB] . 

Let us explain our model and results in more detail. We describe the BEC on a general spacetime metric following 
references [TZHIH]. In the quantum hydrodynamic regime, the BEC phase is described by a mean field classical 
background T plus quantum fluctuations II [20]. The dynamics of is governed by the Gross-Pitaievskii equation. 
At thermal energies much lower than the chemical potential, we can neglect the thermal depletion and assume that all 
the atoms remain in the ground state. Under these conditions, the effect of a gravitational wave on 4/ is just to induce 
a phase shift [2lj. We will discuss this in more detail below. By now, let us focus on the dynamics of the quantum 
fluctuations. For length scales larger than the so-called healing length, n behaves like a phononic quantum field on a 
curved metric. Indeed, as long as the fluctuations in density along the condensate are so small that we can neglect the 
quantum pressure term, the field obeys a massless Klein-Gordon equation nil = 0 where the d’Alembertian operator 
□ = 1 /a /—Bd a {\/ ZI BB ab db) depends on an effective spacetime metric -with determinant g- given by UZHEl 
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(i) 


The effective metric is a function of the real spacetime metric g a ^ -that in general may be curved- and background 
mean field properties of the BEC such as the number density no, the energy density po> the pressure po and the speed 
of sound c s := c^/dp/dp. Here p is the total pressure, p the total density and V a is the so-called 4-velocity flow on 
the BEC, given by the gradient of T. The pressure p and the density p differ from the their bulk counterparts only in 
a small linear perturbation. This description stems from the theory of linearised perturbations of fluids in a general 
relativistic background HU, and thus is valid as long as the BEC can be described as a quantum fluid, that is as long 
as it remains within the quantum hydrodynamic regime m • In the absence of background flows and considering a 
single spatial dimension, we obtain V a = (c, 0). In this case the effective metric reduces to 
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In the absence of a gravitational wave, the real spacetime metric is g a b = g a b^ where = (—c 2 ,1) is the flat two- 
dimensional Minkowski metric and the effective metric Q a b is then a Minkowski-like metric with the speed of light 
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being replaced by the speed of sound c s . The solutions of the Klein-Gordon equation are then massless excitations 
propagating with the speed of sound c s . Therefore, the frequency of the mode uj k is given by the dispersion relation 
ujk = c s |k|, where k is the mode’s momentum. This linear dispersion is valid as long as hk « mo q, where mo is 
the mass of the BEC’s atoms. 

We consider that the BEC is contained in a 1-dimensional cavity trap. Therefore, we choose to impose close to 
hard-wall boundary conditions l22ll24l that give rise to the spectrum, uj n = n7T L Cs ; where L is the cavity length and 
ne {1,2...}. 

The phononic field II(£, x) is then quantised by associating creation and annihilation operators a\ and a k to the 
mode solutions of the Klein-Gordon equation in the effective metric 0 ED], and can be expanded as II(£, x) = 
J2 k [</>&(£, x)a k + </>£(£, x) ajj. The bosonic operators a k and a\. obey the canonical commutation relations. 

By restricting our analysis to phononic Gaussian states, we are able to use the covariance matrix formalism to 
describe the dynamics of II. Gaussian states of bosonic fields and their transformations take a very simple form in the 
covariance matrix formalism. This simplifies the application of quantum metrology techniques to relativistic quantum 
fields [251126] . Considering a collection of N bosonic modes, we define the quadrature operators X 2 n -i = ~^( a n + a n) 

and X 2 n = ~^ 2 i^ a n ~ a n ) wh ere n — 1,...,iV, which correspond to the generalised position and momentum operators 
of the field, respectively. In the covariance matrix formalism Gaussian states are completely defined by the field’s 
first and second moments. Furthermore, quadratic linear unitary operators, such as Bogoliubov transformations, are 
represented by symplectic matrices S that satisfy S T ft S = ft. Here, the matrix ft is the symplectic form defined by 
^ = ©Li n *> n k = —i(Jy and a y is the corresponding Pauli matrix. The first moments of the state are ( X {) and 
the second moments are encoded in the covariance matrix a defined by = (XiXj + XjXi) — 2 (Xi)(Xj). Without 

loss of generality, we restrict our analysis to initial Gaussian states with vanishing first moments, i.e. (JQ) = 0. 

In this work we will consider the following initial state of the field a 0 : 

* 0 = S T isS (3) 

where u = fliag, ;/„, v m ,;/,,,} is called the Williamson form and ig, = coth( ^] k ' T ) are the symplectic eigenvalues of 
the initial state cr 0 (i.e., the eigenvalues of the matrix |if2cro|). We will consider the quantum regime huo k » ksT 
so we can expand u k as follows: 


coth(^) = 1 + 2e~ 2f3k + 0{e~^ k ). 


(4) 


where j3 k := ^ 1 - 

For relatively high frequencies, such as = 2tt x 5 • 10 3 Hz- which corresponds to L = 1/rm and c s = 10 mm/s- 
the above condition entails that we can consider temperatures up to 150 nK -which corresponds to e 2/3 ^ 10. This is 
a relatively high temperature for a BEC where T can be even lower than 1 nK HUES- Note however that we should 
be cautious in extending our analysis beyond a few nK , since we are neglecting thermal depletion in the condensate 
bulk, that is temperature is much smaller than the chemical potential ksT « fi. For typical values of the chemical 
potential ti/ks > 100 nK, which implies that it is reasonable to consider temperatures as large as 10 nK. As the 
temperature grows, the effect of the thermal cloud on the dynamics of the quantum field modes might become relevant 
via the Beliaev damping mechanism [29], which we are not considering here. 

Note that if m > n then, /3 m = m /3 n /n. This means that for m = 2 and n = 1, it is possible to neglect when 
compared to e~^ n . In the regime of temperatures considered here, we obtain v = 1 + j/ 2 ) an q therefore, we can write 
the initial state as 


(To = S T (1 + 21/ 2 )) S, (5) 

where the small contribution takes the form v {2) = diag'{.:E n , x n , x rn , x m } and Xk = e~ /3fc . We will also consider that 
the transformation encoded in S' is a two-mode squeezing transformation. This transformation can be implemented in 
the laboratory, for instance, by artificially modulating the length of the trap [30] . or equivalently by modulating the 
interaction strength BM- However, the degree of squeezing achievable by this method would be limited mainly by 
depletion and losses, so perhaps alternate methods must be considered. Following [5], the frequencies of the squeezed 
modes resonate with the frequency of the gravitational wave, which in turn tunes our detector to a particular frequency 
bandwidth. 

Under the action of the gravitational wave the real spacetime metric g yi/ is transformed from the Minkowski flat 
spacetime metric to a spacetime including a small perturbation i.e. g yv = + h yv (see [4]). In one spatial 

dimension and the traceless-transverse (TT) gauge a perturbation corresponding to a gravitational wave moving in a 




transverse direction can be written as, 
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h - (° 0 ^ 

^ ~ \0 h+(t)J’ 

where h+(t) is typically modelled as a sinusoidal function h+(t) = esin Ut, and e and D are the amplitude and 
frequency of the spacetime ripple, respectively. The change of the real spacetime metric g^ v induces a change of the 
effective metric This in turn generates a Bogoliubov transformation S on the quantum field. The flat- spacetime 
field operators ak are transformed into, = JU [a^aj + ?kj a ]) ? where a\. and d/~ are creation and annihilation 
operators associated to the mode solutions in the perturbed spacetime and akj(h+(t)) and f3kj(h+(t)) are Bogoliubov 
coefficients that depend on the wave’s spacetime parameters. They were computed in [5] in the case of a box-like 
BEC trap under the following assumptions. 

We first consider that beam-pointing laser noise is negligible in the range of kHz frequencies [3TH33] and thus the 
trap can be considered as rigid, and the intensity of the gravitational wave is so small that it remains rigid under 
its action. Second, the frequency of the gravitational wave matches the sum of two modes of interest U = uj m + uj n . 
Third, the BEC-wave interaction time t is long enough uj\ t » 1. The latter two conditions hold for typical sources 
of gravitational waves. Under all the above conditions, the effect of the spacetime ripple is equivalent to a two-mode 
squeezing operation characterised by the coefficient /3 mn , which grows linearly in time in a phenomenon resembling 
the Dynamical Casimir Effect. 

The Bogoliubov transformation generated by the gravitational wave is encoded in the symplectic matrix S e and 
the final state after the transformation is 


where 


S e = 


and the 2x2 matrices Ad mn are given by 
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Here 5ft and ^ denote the real and imaginary parts, respectively. The Bogoliubov coefficients and thus, the final state 
of the field, depend on the amplitude of the gravitational wave e. 

The main goal of this work is to analyse the impact of the temperature T in the bound on the optimal precision 
that can be achieved when estimating the wave amplitude e through measurements on the state cr e . The quantum 
Cramer-Rao theorem states that the error in the estimation of the parameter e is bounded by ((Ae) 2 ) > , where 

H e is the Quantum Fisher Information (QFI) and M the number of probes. The QFI can be computed using the 
Uhlmann fidelity T between the state cr e and a state cr e+rfe with an infinitesimal increment in the parameter. In 

particular one has H e = 8 ( 1 V^C^o-e+dc)) 


de 2 

Now let cr e be a two-mode Gaussian state with zero initial first moments. The fidelity is then given by m 

&e-\-de) 


1 


Va + Vt- + - a’ 


(7) 


where we have introduced T = ^det (id cr e in cr e +de + 12 x 2 ), A = Ardet(yf2cr e + 12 x 2 ) det (ifl cr e+ de + 12 x 2 ) an d 
A = ygdet(cr e + cr e+ d e ). In [26], it is shown that in the case in which the initial state a 0 is pure, V = A + 0(e 4 ) and 
A = 0(e 2 (e + de) 2 ) and so the only relevant contribution to the QFI comes from the computation of 1/VV. In the 
case that we are analysing here, namely a quasi-pure initial state where temperature is only a small perturbation, 
we find that this conclusion is still true. To lowest order, it is possible to show that T = A + 0((x n + x rn )(e + de) 2 ) 
and A = 0((x n + x m ) 2 (e + de) 4 ). This implies that the total contribution of the second square root of 0 cannot 
provide the terms proportional to din? which are the terms that contribute to the final QFI. Furthermore, since 
A = 0((x n +x m ) 2 (e + de) 4 ), one concludes that does not contribute to the final dh 2 term either. Therefore, the 









5 


a) 


Ae 

Vn 


1.5 xlO- 25 


(HzT) 1.X 10- 


5. x 10 -26 



b) 


7 = ( 


c) 


Ae 

Vn 


(HzT) 


500 


1000 


1500 


2000 



500 


1000 


1500 2000 


FIG. 2: Time evolution of the optimal bound for the strain sensitivity as provided by the QFI for n = 1, m — 2 , T = 0 (blue, 
solid), T = 150 nK (red,dashed), m = 6 , T = 0 (black, solid), T = 150nK (green,dotted) and a) r = 10, = 5 10 3 Hz, b) 

r m 10, u>i = 5 10 2 Hz, c) r = 2, wi = 5 10 3 Hz. 


only contribution to the QFI is given by T, as argued before. If one wishes to find at which order the temperature will 
contribute, one needs to compute higher order corrections to all terms in the fidelity. It is possible to show, however, 
that the temperature will contribute linearly to the zero order QFI with a term proportional to x n ,x m . 

We proceed to illustrate numerically what we described above. By using Eqs. 0 and 0 and taking into 
account the discussion below the latter, we compute the QFI H e . Since the analytical computations are cumbersome 
we show here plots obtained by numerical means. In Figj2]we show that the the optimal bound for the strain sensitivity 
provided by the QFI is very robust to initial temperatures. Actually, it is even improved at relatively high temperatures 
-a fact that is related with the already well-known increase of quantum correlations with temperature [351. However, 
note that strictly speaking at those high temperatures the dynamics of the condensate would be different and should 
include damping mechanisms that we are not considering here f36]. The only reason to consider temperatures as 
high as 150 nK is to illustrate the fact that at lower temperatures the plots would be indistinguishable from the 
ones at T = 0. Thus, restricting ourselves to the region ksT « /q where we can neglect thermal depletion in the 
atomic cloud, we conclude that thermal noise in the phononic state does not affect significantly the performance of a 
BEC-phononic gravitational wave detector. A remarkable feature of our mechanism is that the sensitivity improves in 
the kHz regime, as can be seen in Figj3J As discussed in the introduction, this would enable the test of gravitational 
waves coming from binary mergers, which are now beyond the reach of current detectors. Finally, Fig. [2] shows that 
the correction to the vacuum results is, to very good approximation, linear (corresponding to a translation of the 
plot). This corroborates the theoretical analysis which finds that the corrections to the QFI for small temperatures 
corresponds to a small correction to the vacuum QFI. 

Let us now review the validity of some underlying approximations. We have considered that thermal depletion 
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FIG. 3 : Optimal bound for the strain sensitivity as provided by the QFI vs frequency for t = 10, n = 1, ra = 2 , T = 0 (blue, 
solid), T — 150 nK (red,dashed), m = 6 , T = 0 (black, solid), T = 150nK (green,dotted) and a) r = 5, b) r = 7, c) r — 9 . 


is negligible in the atomic bulk and is not significantly modified by the action of the gravitational wave. Indeed, 
in the case of negligible thermal depletion, the wave function of the condensate would acquire a phase T(£) = 
—hk 2 /(2mo) (t — ecos fH/fi), where in a ID box-like potential k = tt c/L and mo is the mass of the BEC’s atoms 
[21] . Since the velocity flows are defined by V^ = cu^/\\u\\, where u^ = h/m^d [18], this means that V = (c, 0) 
both in the absence of the gravitational wave and under its action -in agreement with our assumptions. Moreover, 
since the spacetime ripple only generates a phase shift, it does not change the number of condensed and depleted 
atoms. So if the condensate is initially prepared in a state where thermal depletion is negligible, the thermal depletion 
would remain negligible when the condensate undergoes interaction with the wave. Finally, it is interesting to analyse 
whether the phase shift of the condensate bulk could be used in order to detect the spacetime ripple or not. Indeed 
the QFI H e associated to a quantum state (j>(t) = (j >o is given by H e = \d e ^(t)\ 2 . Therefore, in this case 

H e = (hk 2 /(2mofl) cosOt) 2 . Thus, we first note that, unlike the quadratic growth in time of the QFI in our scheme, 
the QFI of the bulk would only oscillate in time. Furthermore, considering the mass of 87 Rb and the same values of 
L and ft that we are considering in the rest of the work, the maximum value of this atomic QFI is H e ~ 10 —2 , which 
would provide a bound for the strain sensitive of approximately 10 -8 Hz" 1 / 2 with the same parameters of Fig. §i). 
Therefore, the atomic phase shift is extremely less sensitive than our mechanism and cannot be used for gravitational 
wave detection. 

Finally, in order to compare with the sensitivity of state of the art technology (i.e., LIGO, advanced LIGO, VIRGO 
and others), we plot the characteristic strain against the frequency together with that of other detectors. The results 
can be seen in Fig. Q. 

Summarising, we have shown that the performance of our scheme for BEC-phononic gravitational wave detection 
is not significantly affected by the presence of initial thermal noise in the state of the phonons. The gravitational 
quantum resonance between the spacetime ripple and the quasiparticle modes in the BEC trap, is still present in a 
wide regime of temperatures, well within reach of cutting-edge cold-atoms technology. This represents a step further in 
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Frequency / Hz 


FIG. 4: Scaling of the strain sensitivity of our system (solid red line) as a function of frequency in comparison with other 
gravitational wave detectors (solid black lines). 


the feasibility analysis of this novel scheme of gravitational wave astronomy, aiming to complement the open ambitious 
quest for gravitational waves. 


Appendix 

Throughout this work, we have used the QFI as the figure of merit for sensitivity. The QFI provides an ideal 
bound, since it is obtained from an optimization of the classical Fisher information over the set of all the possible 
measurements. However, the optimal measurement might not be easily implementable in the lab and the classical 
Fisher information of accessible measurements could be far from the ideal bound. This could question our choice of 
the QFI as a suitable figure of merit. However, in this Appendix we compute the classical Fisher information of a 
particular realistic measurement and show that it is close to the QFI in a relevant experimental parameter regime. 

We consider that the estimation of the amplitude of the wave is achieved through the measurement of the number 
of phonons in each mode. In [30] . this is achieved for modes carrying different momenta. This scheme relies on 
two key elements: first, a phonon evaporation process, according to which a phonon evolves to an atom of the same 
momentum when the BEC trap is released under adiabatic conditions. Second, after some free-fall time the atoms 
are collected on an array of position-sensitive detectors, from which is possible to infer the arrival time and initial 
velocity -and thus the momenta of the phonons. The frequency range in m is the same that we are considering 
in this work -namely, KHz. 

Then, the classical Fisher information is given by: 

p V" | d-P(e|jO |2 

6 ^P(e\jY de 1 

where P(e\j) is the probability distribution of the state, conditioned to the result j of the measurement- which in this 
case means the detection of j phonons in each mode. In the case of the thermal two-mode squeezed state considered 
in this work, we have: 


( 8 ) 


pm 


1 - e-P" - e-P™ 
cosh (r + /3 nm (e)) 


tanh- 7 (r + /3 nm (e)), 


(9) 


where all the magnitudes have been defined in the main text. Using Eqs. (|8| and § we can compute the classical 
Fisher information and compare it to the QFI. In Fig. <§ we see that the classical Fisher information is actually 
very close to the QFI for relevant experimental parameters. We have observed a similar behavior in all the parameter 
range explored in this work. Therefore, the use of the QFI as a figure of merit is well justified. 
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FIG. 5: Inverse of the classical (blue, solid) and quantum (red, dashed) Fisher information vs time for n = 1, m — 2, 
(jj± — 510 3 Hz, T — 50 nK and r — 2. 
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